Vibrational dynamics and boson peak in a supercooled polydisperse liquid 
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Vibrational density of states (VDOS) in a supercooled polydisperse liquid is computed by diag- 
onalizing the Hessian matrix evaluated at the potential energy minima for systems with different 
values of polydispersity. An increase of polydispersity leads to an increase in relative population 
of the localized high-frequency modes. At low frequencies, the density of states show an excess 
compared to the Debye squared-frequency law, which has been identified with the boson peak. 
The height of the boson peak increases with polydispersity. The values of the participation ratio 
as well as the level spacing statistics demonstrate that the modes comprising the boson peak are 
largely delocalized. Interestingly, the intensity of the boson peak shows a rather narrow sensitivity 
to changes in temperature and is seen to persist even at high temperatures. Study of the difference 
spectrum at two different polydispersity reveals that the increase in the height of boson peak is 
due to a population shift from modes with frequencies above the maximum in the VDOS to that 
below the maximum, indicating an increase in the fraction of the unstable modes in the system. 
The latter is further supported by the facilitation of the observed dynamics by polydispersity. Since 
the strength of the liquid increases with polydispersity, the present result provides an evidence that 
the intensity of boson peak correlates positively with the strength of the liquid, as observed earlier 
in many experimental systems. 

PACS numbers: 63.50.Lm, 64.70.pm, 82.70.Dd 
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I. INTRODUCTION 



A glass behaves mechanically like a solid but struc- 
turally like a liquid with relaxation times varying from a 
few minutes to several centuries. The elasticity of a solid 
is described in terms of phonons, which are quantized 
vibrational excitations. Propagating acoustic phonon- 
like excitations have been observed in glasses and glass- 
forming liquids down to wavelengths comparable to inter- 
particle distance However, their description is ren- 
dered difficult since glasses lack the translational in- 
variance of crystalline solids. A ubiquitous feature in 
the physics of glasses is the anomalous behavior of the 
low-frequency part of the vibration spectrum and the 
corresponding thermal properties While the ori- 

gin of the linear low-temperature specific heat is com- 
monly attributed to the existence of double-well poten- 
tials or two-level systems, there is a considerable debate 
about the so-called boson peak [1, [1]. This peak shows 
up in the vibrational density of states(VDOS), g{i') as 
an excess contribution, compared to the usual Debye 
behavior[(7(i^) oc i/'^]. It is called boson peak because the 
temperature dependence of its intensity scales roughly 
with the Bose-Einstein distribution. In addition to the 
presence of boson peak (BP) (or excess density of states) 
at low frequencies, one also observes a high-frequency 
exponential tail [s] in the reduced density of states spec- 
trum. 

The interpretation of the boson peak has been a chal- 
lenge to experimentalists and theoreticians and is a sub- 



ject of controversial discussions. While some authors at 
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tribute it to local/quasi- local vibrations 
others attribute it to collective motions [5, 
mechanism of the BP formation in glasses was proposed 
based on the concept of interacting quasi-local oscilla- 
tors Q. Boson peak has also been explained in terms 
of the afline-non-affine crossover at a certain mesoscopic 
length scaleflO]. It has also been interpreted as the signa- 
ture of a phase transition in the space of the stationary 
points of the energy, from a minima-dominated phase 
(phonons) at low energy to a saddle dominated phase 
(without phonons) 11]. The boson peak has also been 
linked to those motions giving rise to the two-level-like 
excitations seen at still lower temperatures [12]. Recently 
evidence was presented from numerical studies suggestive 
of the equality of the boson peak frequency to the loffe- 
Regel limit for transverse phonons above which trans- 
verse phonons do not propagate, and the boson peak 
was attributed to transverse vibrational modes associ- 
ated with defective soft structures in the disordered state 

Optical heterodyne-detected optical Kerr effect data 
on supercooled acetylsalicylic acid and dibutylpthalate 
display highly damped oscillations with a period of a 
few pico seconds as the temperature is reduced to and 
below the Mode Coupling Theory temperature, Tmct 
p^ . The authors interpret this as the time domain sig- 
nature of the boson peak and explain that the increased 
translational-rotational coupling is responsible for the 
boson peak as T < Tmct- The main experimental tools 
used for the investigation of VDOS and boson peak in 
glasses are Raman scattering and the quasi-elastic neu- 
tron scattering(QENS). Recently, single-molecule spec- 
troscopic studies have also emerged as a potential tool to 
study the boson peak in glasses [3]- The main difficulty 



here is the conflicting interpretations regarding the origin 
of the boson peak. For instance, nuclear inelastic scat- 
tering studies of ferrocene probe molecules in toluene, 
ethylbenzene, dibutylphthalate and glycerol glasses show 
that a significant part of the modes constituting the bo- 
son peak is of collective character Q. However, from 
the single molecule spectroscopic studies of many single 
tetra-tert-butylterrylene molecules embedded to amor- 
phous polyisobutylene the authors conclude that the low- 
energy vibrational excitations have a local character Q. 

The theoretical as well as experimental studies are 
done on systems/samples in the glassy state. However, 
since the the phonon-like excitations also persist at el- 
evated temperatures in the supercooled liquid regime 
[l|, once can study the vibrational dynamics of super- 
cooled liquids via molecular dynamics simulations. Vari- 
ous groups have studied the vibrational dynamics of su- 
percooled liquids [H, [IB, [13, [ll], where the focus has 
been on the characteristics of the high-frequency or low- 
frequency vibrational modes. The main emphasis of 
these studies have been on the collective/local nature of 
the vibrational modes. 

In this study we investigate the vibrational dynam- 
ics and boson peak in a polydisperse Lennard- Jones liq- 
uid. Polydisperse liquid is one of the simplest model 
systems that exhibit glass transition and can be conve- 
niently studied via both experiments [Tol . [20| and com- 
puter simulations as the size distribution of particles pre- 
vents crystallization [21, 22,|23|. The rest of the paper is 
organized as follows. In section |TT] we describe the model 
and computational details. In section Hill we present our 
results and give detailed discussions on the same. We 
give our concluding remarks in section [TVl 



II. COMPUTATIONAL DETAILS 

Micro canonical ensemble molecular dynamics (MD) 
simulations are carried out in three dimensions on a sys- 
tem of Lennard- Jones (LJ) particles of mean diameter 
with polydispersity in both size and mass. The polydis- 
persity in size is introduced by random sampling from 
the Gaussian distribution of particle diameters cr, 



1 , 1 cr — (T,2i 



(1) 



The standard deviation S of the distribution divided by 
its mean W gives a dimensionless parameter, the polydis- 
persity index, S 



fixed volume fraction, (/) — 0.54. Three different system 
sizes were chosen, N = 256, 500 and 864. The results are 
found to be qualitatively the same for the three different 
system sizes studied. 

The interactions between the particles are given by the 
shifted-force LJ potential 



where i and j represent any two particles and 



(4) 



(5) 



The LJ interaction parameter is assumed to be the 
same for all particle pairs and set equal to unity. The 
particles are enclosed in a cubic box and periodic bound- 
ary conditions are used. The cutoff radius is chosen 
to be 2.5ct. The time step used for integrating the equa- 
tions of motion is 0.001. All quantities in this study are 
given in reduced units (length in units of cr, temperature 
in units of and time in units of r = 

The vibrational density of states is obtained from the 
Normal Mode Analysis by solving the secular equation, 



^ 



(6) 



Here F is the mass-weighted potential energy matrix 
(also known as the Hessian matrix), 



F ■ — 



(7) 



Equationinican be solved to yield a set of eigen values (the 
square of the vibrational frequencies, v'^) and correspond- 
ing eigen vectors (normal mode displacement vectors), e^. 
From the equilibrium liquid configurations generated by 
the MD simulations, one constructs the potential energy 
minima or the inherent structure via conjugate gradi- 
ent minimization. From the potential energy minima, 
one generates the Hessian matrix, the diagonalization of 
which would yield the eigen values and the eigen vectors. 
The normal modes thus obtained are called the quenched 
normal modes (QNM). 

Instantaneous normal modes (INM) can be obtained 
by diagonalizing the Hessian constructed for the equi- 
librium liquid configurations. Since the instantaneous 
liquid configuration is not necessarily a potential energy 
minimum, one gets both unstable modes (negative eigen 
values or imaginary frequencies) and stable modes (pos- 
itive eigen values). 



S ■ 



The mass of particle i is scaled by its diameter. 
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III. RESULTS AND DISCUSSION 



A. Polydispersity Effects on Vibrational Density of 
States 



We have chosen m — 1.0. The simulations are carried out 
at different values of the polydispersity index, S but at 



The configuration-averaged vibrational density of 
states, g{iy) for the quenched normal modes is shown in 
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FIG. 1: (a)Quenched normal mode spectra at different values 
of r. The data is for 5" = 0.20 and system size N = 256. 
(b)Instantaneous normal mode spectra for different values of 
T. As T decreases the INM spectra approaches closer to the 
quenched spectrum (dark circles). 
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FIG. 2: (a) g(u) for S = 0.10 and S = 0.20 systems and 
their difference (inset) (b) Same as (a) but for S = 0.10 and 
S = 0.20 systems without mass polydispersity. Data shown 
for a system of = 256 particles. 



FIGURE ma) for S = 0.20 at different values of tem- 
perature. The width of the bin chosen to build the his- 
togram is 0.50. The VDOS is rather featureless as has 
been observed by Rahman et al for single-component 
Lennard- Jones system. From FIGURE [Ija), we see that 
the quenched normal density of states are only weakly 



sensitive to the temperature of the parent liquid. The 
rather narrow sensitivity of the quenched normal mode 
spectra to temperature is to be contrasted with the in- 
stantaneous normal mode spectra obtained from equi- 
librium liquid configurations (See FIGURE [T]Jb)). The 
imaginary modes in the INM spectra are displayed on 
the negative side of the frequency axis. Both the stable 
(positive frequencies) and unstable (imaginary frequen- 
cies) modes of the INM spectrum show a pronounced 
change with temperature. As the temperature decreases, 
the fraction of the unstable modes decreases whereas 
that of the stable modes increases. At low temperature 
(T < 0.40), the fraction of the unstable modes is con- 
siderably less, indicating that most of the particles are 
located near the potential energy minima most of the 
time. 

Polydispersity has three noticeable effects on the vi- 
brational density of states. Firstly, as polydispersity in- 
creases, the number of low frequency modes increases 
whereas the number of high frequency modes decreases 
(See FIGURE[2|). Because of these compensating changes 
occurring in the high and low frequency regions, one 
observes a crossover in the density of states between 
S = 0.10 and S = 0.20 systems. For the data shown 
in FIGURE [2l this happens at a frequency ly ~ 12.0 for 
systems with mass polydispersity and v ~ 16.0 for sys- 
tems without mass polydispersity. 

Secondly, when mass polydispersity is present, there is 
a second crossover point in the density of states between 

= 0.10 and S = 0.20 systems at a frequency, v ^ ^2 
(FIGURE [21a)). For frequencies higher than this value, 
there is an excess of high-frequency modes that increases 
with S. This is best seen in the semi-log plot of g{v) in 
FIGURE [21 The plot clearly shows that as polydispersity 
increases there is a substantial increase in the number of 
high frequency modes. The polydispcrsity-dependence of 
the INM spectra is shown in FIGURE [H The S = 0.20 
system has relatively higher number of high frequency 
modes for both the stable and unstable branches. In FIG- 
URE [5l we show the semi-log plot of VDOS for different 
S but without mass polydispersity (i.e. all masses set 
equal to unity, — 1.0). As we can see from the figure, 
the number of high-frequency modes decreases with S 
for systems having no mass polydispersity. Thus the ex- 
cess high frequency modes (whose fraction increases with 
increase in polydispersity) is due to the mass polydisper- 
sity effect rather than size polydispersity effect. (Below 
we show that the high frequency vibrations are all local- 
ized.) Thirdly, the vibrational density of states spectrum 
becomes narrower with polydispersity. In FIGURE [Bj we 
plot the full width at half maximum (FWHM) of g{v) for 
S = 0.10 and 5* = 0.20 as function of T. The FWHM 
decreases sharply with S but shows only a weak increase 
with temperature. 

The localization properties of the normal modes can 
be quantified via the participation ratio (PR) which is 
a measure of the number of particles participating in a 
given vibrational mode. The participation ratio of mode 
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FIG. 3: Semilog plot of density of states. The data is for 
T = 0.50. The high frequency modes are localized (See Fig 



FIG. 5: Semilog plot of density of states for different S with- 
out mass polydispersity. The data is for T — 0.50. 
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FIG. 4: Instantaneous normal mode spectrum for 5 = 0.10 
and S = 0.20 systems at T = 0.80. 



is defined as 



PR^ 



3N 
a=l 



(8) 



The participation ratio is of the order of jj for locaHzed 
modes and is of the order of 1 for extended modes. In 
FIGURE [3 the participation ratio of the modes is plot- 
ted. The averaging is done over modes corresponding 
to eigenvalues in a histogram bin of width 0.50. As one 
can see from the plot, for v > 30.0, the participation ra- 
tio is very low, indicating that the high-frequency modes 
are all localized. Therefore the high-frequency tail of the 
normal mode spectrum can be attributed to the localized 
vibrations. In a crystalline solid, localized vibrations oc- 
cur due to the presence of very light impurity atoms or 
interstitial atoms (that cause large lattice strains) . It has 
been shown that when the mass of the impurity atom in 
a crystal is smaller than that of the other particles, the 
vibrations do not propagate through the system but gets 
localized around the impurity particle [25| . 

In the frequency range 3.0 < i' < 35, one finds delo- 
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FIG. 6: Full width at half maximum, FWHM for S = 0.10 
and 5* = 0.20 as a function of T. 



calized modes as implied by the high values of PR. The 
participation ratio of the modes decreases with polydis- 
persity for modes of all frequencies (See FIGURE [7|) im- 
plying that the number of particles participating in a 
mode of given frequency decreases with S. Thus the vi- 
brations become more localized with S. This means that 
as size disparity among the particles increases, the system 
cannot sustain propagating modes. It is interesting to 
note here that the large size disparity among the particles 
also leads to the suppression of growth of dynamic het- 
erogeneity with polydispersity in supercooled liquids(See 
[2^). Dynamic heterogeneity in supercooled liquids in 
its simplest sense means clusters of fast-moving particles 
that move together for a certain amount of time before 
they get decorrelated. When size disparity is large, the 
particle motion gets decorrelated much faster and hence 
the formation of dynamic clusters is suppressed at higher 
polydispersity. In the inherent structure formalism devel- 
oped by Stillinger and Weber [2^, [l^ , the configuration 
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FIG. 7: Participation ratio (PR) of the modes for system 
sizes, (a) A'^ — 256 and (b) A'^ = 864. PR gives the number of 
particles participating in a given mode. The data is shown for 
different S at T = 0.50. The plot shows that a S increases, the 
number of particles participating in a given mode decreases 
for all frequencies. 
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FIG. 9: Excess density of states (Boson Peak), versus ly. 
Data shown for different S with the temperature of the parent 
liquid at r = 0.50 and for system sizes, (a) N = 256 (b) N 
= 256 but without mass polydispersity (c) N = 500 and (d) 
N=864. 
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FIG. 8: Basin curvature as a function of temperature for S - 
0.10 and S = 0.20 systems. 




FIG. 10: Temperature dependence of the Boson Peak. The 
data is shown for S = 0.10 for different temperatures of the 
parent liquid. 



space of the liquid is divided into basins of local potential 
energy minima(or inherent structures). In the harmonic 
approximation, which is valid at sufficiently low temper- 
atures, each basin is treated as a harmonic well. In FIG- 
URE. [8] we plot the quantity N^^ X]fc=i ^og{hvk) as a 
function of T for S* = 0.10 and S = 0.20 systems. This 
quantity is an indicator of the average curvature of the 
basins |28j . Since it is a sum of logarithms it is very sensi- 
tive to the tail of the spectrum. From the plot we see that 
the average curvature of the basins increases with T but 
decreases sharply with polydispersity. In other words, in- 
creasing polydispersity leads to a flattening of the basins 
which would in turn facilitate inter-basin transitions and 
thus enhances the mobility of the s yste m. The latter has 
been termed the lubrication effect [23| . 



B. Boson Peak 

We plot the reduced density of states, g{v)/v'^ in FIG- 
URE [i for three values of S, namely S = 0.10, 0.15 and 
0.20. All the three systems exhibit the boson peak fea- 
ture viz. excess density of states as compared to the 
prediction of the Debye model {g{v) ^ v^). As seen from 
the figure, the intensity of the boson peak, Ibp increases 
with S for all the three system sizes studied. However, 
there is no noticeable change in the frequency of the Bo- 
son peak, vbp- The boson peak feature is seen even when 
we switch off the mass polydispersity (See FIGURE [9] 
(b)) implying that the size polydispersity effect alone can 
give rise to the observed features. Furthermore, the bo- 
son peak feature is seen even for the quenched normal 
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FIG. 11: The level spacing distribution, P(s) for S = 0.10 
and S = 0.20 systems for the full eigen value spectrum. Data 
is shown for 5" = 0.10 and S = 0.20 systems at T = 0.50. 
The symbols are the data points and the thick lines are fit to 
the data according to the equation, P{s) = As^exp{—Bs^). 
The values of f3 are 1.02 and 1.08 for S = 0.10 and S = 0.20 
systems, respectively. 



modes obtained from equilibrium configurations at very 
high temperature (i.e. T ^ Tmct), as shovifn in FIG- 
URE [TOl As T increases, the boson peak height increases 
but the peak shifts to lower frequencies. Such a trend has 
been observed in experiments as well [1^, [s^l- However, 
the temperature-dependence is weak here. The boson 
peak feature at high temperatures shown in FIGURE 
[TUl is to be contrasted with the observation made in an 
earlier paper (ll| where the authors studied the appear- 
ance of the boson peak in soft sphere binary mixture. In 
their study, the authors found that as T increases above 
Tmct, the boson peak feature disappears. The authors 
interpret the boson peak as a manifestation of an under- 
lying crossover of the parent liquid's configuration from 
a saddle dominated dynamics to a minima dominated be- 
havior. However, the observation of boson peak at high 
temperatures makes such an interpretation questionable. 

The fragility of the current model polydisperse super- 
cooled liquid decreases with polydispersity index, S ^2Z| . 
This suggests an inverse correlation between the fragility 
of the liquid and the intensity of the boson peak. Such 
a correlation has actually been observed previously for 
many liquids in experiments [3l| as well as in the simu- 
lation studies of model glass formers [l^. Shintani and 
Tanaka carried out extensive simulations on 2D spin 
liquid glass former in which the fragility can be varied 
over an extremely wide range. This was achieved by 
varying the anisotropic part of the potential, A which 
is a measure of the strength of frustration against crys- 
tallization. As already mentioned, the authors found an 
inverse correlation between the fragility and the inten- 
sity of the boson peak. In the current model system the 



polydispersity index, S plays the same role as A. Auer 
and Frenkel have shown that crystal nucleation in a poly- 
disperse colloid is suppressed due to the increase of the 
surface free energy ^] and studies by several groups [s^l 
show that the glassy amorphous phase becomes the equi- 
librium phase beyond a terminal value of polydispersity. 

If a correlation between the fragility and the intensity 
of the boson peak exists, then it would mean that there 
is also a link between the fast intra-basin vibrational dy- 
namics and the slow inter-basin diffusive dynamics (See 
[3^). However, it appears hard to reconcile such a link 
with the Adam-Gibbs paradigm as the latter predicts 
a relationship between the relaxation time, r and the 
configurational entropy, Sc (t ~ ea;p(-j^)). Clearly, on 
the time scale typical of vibrational dynamics (few pico 
seconds), the system would not have sampled sufficient 
configurations as would be necessary to define Sc- Hence 
it is difficult to understand the correlation between the 
intensity of the boson peak and fragility from this per- 
spective. 

It is interesting to note in this connection that Angell 
et al [3^ have already suggested that the Boson peak 
can serve as a signature of configurational excitations of 
the ideal glass structure i.e. the topologically diverse de- 
fects of the glassy solid state. This means that the boson 
peak is related to the topographical features of the poten- 
tial energy landscape and is thus involved in determining 
the fragility of the liquid. As already discussed in Sec- 
tion [TTTXl at higher polydispersity we have flatter basins 
(FIGURE [5]). This is consistent with the decrease in the 
fragility with S, as fragile liquids have rugged heteroge- 
neous landscape whereas strong liquids have smoother 
landscapes, in accord with their constant activation en- 
ergy predicted by their Arrhenius behavior [36j . 

The strong liquids in Angell's fragile/strong classifica- 
tion are usually network glass formers like Si02, Ge02 
etc., and they appear to lie almost on the opposite spec- 
trum of our polydisperse liquid system. For a polydis- 
perse LJ liquid, the decrease in fragility with polydis- 
persity is via dynamic facilitation by smaller particles 
[22I |23[ . An obvious manifestation of the dynamic facili- 
tation is the size-dependent glass transition temperature 
|.22j ; as temperature is lowered the larger particles freeze 
in first, followed by the smaller ones. Therefore dynamic 
facilitation by polydispersity implies that at higher poly- 
dispersity not only the system has smaller barriers to dif- 
fusion but it also has more relaxation channels available 
to it. 

Here it is interesting to note that Shirmacher at al Q 
have shown that if a system of coupled harmonic oscil- 
lators (with spatially fluctuating nearest-neighbor force 
constants on a simple cubic lattice) is near the borderline 
of stability a low-frequency peak appears in g{y)lv^ as a 
precursor of the instability. In their model system, when 
the amount of the negative force constants becomes too 
large, the system becomes unstable and the boson peak 
feature shows up. Furthermore, as the fraction of nega- 
tive force constants increases the peak intensity increases 




the case of delocalized states, according to the Gaussian 
orthogonal random matrix ensemble, we get 



P(s) = Asf^expi-Bs^ 



(9) 



FIG. 12: The level spacing distribution, P(s) for the eigen 
values in the range 3.0 < v < 4.0 (near the boson peak) 
shown in circles and for eigen values in the high frequency 
range (v > 35.0) shown by squares. Data is shown for S = 
0.10 at T = 0.50. The symbols are the data points and the 
thick lines are fit to the data according to the the Gaussian 
orthogonal ensemble (thick line) and Poissonian distribution 
(dashed line) 



and the peak shifts towards lower frequencies. Instanta- 
neous normal mode analysis shows (see again FIGURE[4|) 
that the fraction of unstable modes increases with poly- 
dispersity which implies more pathways for diffusion [l6| 
at higher polydispersity. Thus the increase in the in- 
tensity of boson peak with polydispersity can be under- 
stood in terms of dynamic facilitation by polydispersity. 
The relationship with the strength/fragility of the sys- 
tem, however, still remains a bit unclear. 

One of the key issues in the interpretation of the bo- 
son peak feature is whether the modes comprising it are 
localized or extended. Figure [7| shows that the modes 
at the frequency range where boson peak appears are 
largely delocalized with high values of participation ratio 
[PR > 0.6). For frequencies less than z/^p, the partici- 
pation ratio drops suddenly. This sudden drop of PR at 
low frequencies has been attributed to finite size effects 
p^ . These low- frequency modes have been shown to be 
extended, notwithstanding their low PR values. How- 
ever, note that the frequency of the boson peak coincides 
with the frequency below which PR suddenly drops. The 
reason for this behavior is not clear to us. Furthermore, 
at the boson peak frequency, the PR decreases with S 
indicative of a correlation between localization of modes 
and boson peak intensity. These features are common to 
all the three system sizes (N=256, 500 and 864) studied. 

An alternate way to check whether the boson peak 
is indeed associated with localized modes (or not) is by 
means of level distance statistics d, 113] ■ The level spac- 
ing distribution, P{s) for the random-matrix models is 
defined as the probability of finding the next-nearest- 
neighbor eigenvalue of the spectrum to be at a distance 
s, i.e. s = -^-^ — where D is the mean level spacing. In 



For localized states, one gets a Poissonian distribution, 
P{s) = exp{—s). This is due to the fact that the delocal- 
ized states show level repulsion, whereas localized ones 
do not. In FIGURE [TT] we have plotted the level spacing 
distribution P{s) for the full eigenvalue spectrum along 
with the fit for the orthogonal random matrix model. 
The values of f3 obtained from the fit are 1.02 and 1.08 
for S — 0.10 and S — 0.20 systems, respectively. This 
shows that these systems display the behavior accord- 
ing to the Wigner surmise. In FIGURE [121 we plot the 
level spacing distribution for frequencies near the boson 
peak (3.0 < ly < 4.0) and also for the high frequencies 
{v > 35.0) where the modes are truly localized (See[71). 
From the plot it is clear that statistics at high frequen- 
cies follow Poissonian distribution as expected for local- 
ized states. In the vicinity of the boson peak, however, 
we find a Gaussian distribution according to the Gaus- 
sian orthogonal random matrix ensemble, which means 
that the corresponding states are delocalized. However, 
it should be noted here that the transverse phonons are 
shown to be localized at the bosoiipeak frequency ac- 
cording to the loffe-Regel criterion [l^ . 



IV. CONCLUDING REMARKS 

Let us first summarize the main results of this paper. 
We have computed the vibrational density of states, g{iy) 
in a polydisperse Lennard-Jones system, for three system 
sizes and three different polydispersity. Polydispersity 
is shown to have a significant effect on the vibrational 
density of states. Increase in the polydispersity leads 
to an increase in the localized high frequency modes. 
As polydispersity increases, there is a softening of vi- 
brational modes manifested in the population shift from 
modes with frequencies above the maximum in g[i') to 
that below the maximum. The reduced density of states, 
g{v)/v'^ exhibits the boson peak feature. The boson peak 
is seen to exist even for high temperature liquid config- 
urations and shows a rather narrow sensitivity to tem- 
perature. The intensity of the boson peak increases with 
polydispersity, which can be understood in terms of dy- 
namic facilitation by polydispersity by increase in frac- 
tion of the unstable modes. 

The results clearly show a correlation between an in- 
crease in the peak height of the BP with an increase in 
the strength of the liquid - both are consequences of in- 
crease in polydispersity. This is in agreement with known 
experimental results. We find that while the modes com- 
prising the BP are largely delocalized, there is a sharp 
drop in the PR of the modes that exist just below the 
BP. The last observation could be due to the finite size 
of the system. The delocalized nature of the BP appears 



to be quite robust, within the three system sizes consid- 
ered here, as supported both by the participation ratio 
and level spacing statistics. 

The observed weak temperature dependence of the bo- 
son peak seem to suggest absence of any temperature 
mediated phase transition in the stationary points of the 
energy landscape from a saddle dominated to a minima 
dominated regime. However, all the evidences do sug- 
gest a close relation between boson peak and unstable 
modes of the system. The precise nature of these unsta- 



ble modes, however, remains unclear at this point. 
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